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'sj- , Abstract 

\o , 

■ The mass spectrum is calculated for those dileptons which are produced in the 

. early phase of a heavy ion collisions via the direct production AW — > X and via 

Q>^ \ the Compton process GN — > X with prompt gluons radiated in preceding AW 

' interactions. Both mechanisms produce a mass spectrum which decreases steeply with 

• invariant mass of the pair and which is below the CERES data for Pb-Au collisions 
I , by about one order of magnitude. 

F; ; 1 Introduction 

> . 
•I— I . 

^ . Dileptons with invariant mass below 1 GeV have been measured in proton-nucleus [Q] and 
^ ! heavy ion collisions ^ at the CERN SpS accelerator. For heavy ions (laboratory energy 
158 A GeV) a sizable enhancement of the order of a factor three has been observed over 
what is calculated from the e+e~ decay of the hadrons in the final state. Also the observed 
shape is very different from the calculated one. The data have not yet found an unambiguous 
explanation. Most theories locate the origin for the observed dileptons in the hot and dense 
phase of hadrons 0]-[|T^. This phase is rather late in the time evolution of a heavy ion 
collision. 

In this note, we take the opposite point of view and investigate the mass spectrum of the 
dileptons which arise from the very early stages of the heavy ion reaction in which partons 
are the relevant degrees of freedom. We calculate (i) the direct production AW — ^ l^l^X 
in the light cone approach and in a parton model and (ii) the lepton production via a 
gluonic Compton process GN — > l^l~X from prompt gluons which are radiated in preceding 
nucleon-nucleon interactions. These prompt gluons have been recently identified as as an 
important source for charmonium suppression in heavy ion collisions |jl3|. This observation 
has triggered the present investigation. 
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2 Direct production of dileptons: light cone approach 
versus parton model 

First we discuss the direct production of lepton pairs, which is aheady calculated using a 
parton model and perturbative QCD [Q. However, in the experiments under consideration, 
the invariant masses of the observed pairs are below 1 GeV. In this domain pQCD may be 
questionable and we use a phenomenology based on the light cone approach which includes 
nonperturbative effects [|15|, [16], O. 



(a) 





Figure 1: Diagrams for the two mechanisms for production considered in this paper, (a) 
Dilepton production in NN collisions which can be viewed as bremsstrahlung in qq scattering 
(light cone approach) or as gluonic Compton scattering [Gq '~f*q) from the gluon cloud of 
the second nucleon (parton approach), (b) Dileptons from prompt gluons Gq —>■ '~f*q, where 
the on-mass-shell gluon G (crossed by a dashed line) is produced in another NN collision. 

In the light cone approach, the Drell-Yan type process illustrated by one of the pos- 
sible Feynman diagrams in Fig. |l|a is viewed in the target rest frame, where it looks like 
bremsstrahlung from a NN collision. The incident quark from A^i scatters off a nucleon N2 
in the target and radiates a massive photon 7* which then decays into a lepton pair. The 
cross section takes the form [|15], [19], |l^ 



rfVT|^^.g(a,fY)|Vgg-(arT,s) 



d{lna) 



where a is the fraction of the light cone momentum of the quark carried away by the photon 
and tt is the transverse separation between the quark and the photon. \E'^*q(a,fr) is the 



light cone wave function for the transition q 
transverse (T) and longitudinal (L) photons 



g7*. We give the explicit expressions for 



l^:;*,(a,^T) 
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with 



a 



where rriq is the quark mass and 



Pr 



Note that 



although there is only a single quark in the initial state the cross section (|I|) for radiation 
of a photon depends on the cross section aqq{arT, s) for the scattering of a qq dipole off a 
proton. Here, ary is the transverse separation of the dipole and s is the cm. energy of 



the quark nucleon system squared. This can be understood as follows |T5| , |I6| . The physical 
quark is represented as a coherent superposition of different Fock states. In our case, we take 
only the bare quark q and the g7* states into account. If these two Fock states would scatter 
with the same amplitude off the target, coherence would be undisturbed and no radiation 
produced. The radiation amplitude is proportional to the difference between the scattering 
amplitudes of the two Fock components. The dipole cross section enters, because of the 
interference between the two graphs in fig. (la). The impact parameter of the projectile 
quark serves as center of gravity for the 57* fluctuation. The difference in impact parameter 
between the parent quark and the quark in the fluctuation is ar-r- Therefore, the two graphs 
have a relative phase factor exp(2afY ■ Pt)- The antiquark appears, when one takes the 
complex conjugate of one of the graphs. It is now easy to see, that in the absolute square 
of the two graphs, the color screening factor [1 — exp{iafT ■ pr)] of crqq{c(i^T, s) emerges from 
the four different possible attachements of the two gluons. Therefore the dipole cross section 
appears in (Q), although there is no physical dipole in the system. 

Perturbative QCD predicts aqq{rT, s) oc for — >■ 0. Indeed, a colorless qq dipole can 
interact only via its color dipole moment. However, since we are interested in the low mass 
region, which corresponds to rather large separations r^, we rely on phenomenology and 
employ the modification Ijl^ of the saturation model of 



1 — exp 




(5) 



where tq (s) = 0.88(s/so)~°'^^ fm, sq = 1000 GeV^. This cross section is proportional to 
for rrp —>■ 0, but flattens off at large rj^. The energy dependence correlates with r^. At small 
tt the dipole cross section rises faster with energy than at large separations: 



o-o(s) 



a 



■n-p, 
tot\ 



1 + 



3rl 



(6) 



where cr^^(s) = 23.6(s/so)°'°^ mb and (r'^hln — 0-^4 fm^. With this choice, value and energy 
dependence of the total cross section for pion proton scattering is automatically reproduced. 
This cross section also allows to describe well the structure function F2 in DIS in a wide 
interval of energies and [|17 



The partonic cross section (|l|) has to be embedded into the hadronic process. The cross 
section of direct dilepton production takes the form 



dMdy 



9(y 

3^ Y""" 



1 Xi 

da—^ 



c/o-(gA^ YqN) 
d{\n.a.) 



+ {y^ -y), (7) 



where y is the rapidity of the lepton pair in the cm. frame and Xi = {^Xp + 4M2/s + xf)/2. 
The first term corresponds to radiation from the projectile quarks and the second term to 
radiation from the target quarks (see Fig |I]a), which corresponds to replacement y — —y in 
the first term. Since the dominant contribution to the cross section comes from large values 



3 



of xi/a, we neglect antiquarks in the projectile, fg (xi/a) = 0. We parameterize the valence 
quark distribution in the well known form 

fu{x) = (8) 

Ux) = ^il-x)\ (9) 



where Cu,d are defined by the normalization to the number of valence quarks Nu^d in a nucleon 

1 

dxfu,d{x) = Nu,d. (10) 



In the case of the proton-nucleus {pA) and nucleus-nucleus {AB) collisions the cross section 
of the direct production is given by the integration of the nuclear densities pa,b over the 
impact parameter 

dA^y = JdhJdsljz^pAs,ZA)ljzBPB{h-s,z,) ^ = AB^. (11) 

In this expression we neglect small shadowing effects. To obtain the rate of the dilepton 
events per interaction one has to divide expression (pA]) by the total cross section, which can 
be calculated using the Glauber expressions 

apA = 1 c/b [l - exp (-or]^^TA(b))] , (12) 
aAB = Jdh 1- exp (^-a™^ J rfsT^(s) Tsih - s)) , (13) 

where cr]^^ ~ 30 mb denotes the inelastic AW cross section and we used the standard Woods- 
Saxon parameterization |TH] for the nuclear density pa,b in 

/oo 
dzpAA^^^)- (14) 
-oo 

There are still two sources of uncertainty: 

• The dipole cross section aqg^rx, s) cannot be calculated perturbatively at large trans- 
verse separations, which are important in the small mass region, however. We employ 
the phenomenological expression (^. 

• The light cone wave function \E'^5(a,fY) depends on the mass nig of the quark. To 
illustrate its influence we calculate the cross sections for the two extreme cases of the 
constituent quark mass rrig = 150 MeV and rrig = 300 MeV. 

Our results for proton (p-Be, p-Au) and heavy ion (Pb-Au) scattering are shown in 
Figs. under the label "light cone" . The border lines of the shadowed areas correspond to 
different assumptions about the constituent mass rrig (150 MeV and 300 MeV for the upper 
and lower lines, respectively). We believe, that also a different dipole cross section cannot 
produce values, which are significantly out of the shaded area, because agg is well constrained 
to describe the structure function F2 and hadronic cross sections. 
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Figure 2: Data on dilepton events in p-Be and p-Au scattering together with predictions 
from different mechanisms. The various curves correspond to meson decay channels. The 
shaded areas represent our resuhs for direct production (in the hght cone and the parton 
model approaches) and via prompt gluons (see section |^): the upper limits correspond to 
the quark mass rriq = 150 MeV and the lower ones to rriq = 300 MeV. 

It is instructive to compare these results with predictions from a parton model which 
treats direct radiation of the lepton pair as Compton scattering of a target gluon on a beam 
quark (or vice versa). This way of description should be identical to one employed above 
(as one can see from Fig. |^). However, all phenomenological modifications introduced via 
the effective qq cross sections are missing. The subprocess qG — *• l^l~X is then calculated 
perturbatively with the help of the valence quark and gluon distributions in the proton. 

The main uncertainties in this case originate from: 

• Possible failure of pQCD calculations for the soft Compton amplitude. 

• Poor knowledge of the gluon distribution function, especially at a soft scale. 
We use the following gluon distribution function, 

/G(x) = %(l-x)^ (15) 

where Cq is defined from the condition that gluons carry a half of the nucleon momentum 

/•I 1 

/ dxxfcix) = -. (16) 
Jo 2, 
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Figure 3: Data 0, |[ on dilep- 
ton events in Pb-Au scattering. 
The different curves correspond 
to meson decay channels. The 
shaded areas are our results for 
direct production in the light 
cone approach and via prompt 
gluons (see section H): the upper 
limits correspond to the quark 
mass mq = 150 MeV and the 
lower ones to niq = 300 MeV. 
Both processes contribute sepa- 
rately to the dilepton production 
and therefore their contributions 
have to be added. 



M (GeV/c^) 

The pQCD calculation of the Compton amplitude is similar to that for an interaction with 
on mass shell gluons presented in the next section. 

The results for proton (p-Be, p-Au) and heavy ion (Pb-Au) scattering are shown in 
Figs. ^ 1^. At low M the results of pQCD (labeled as "parton model") are higher than 
those from the light cone calculations but start to agree for M>1 GeV. The reason for this 
is that the saturation of the dipole cross section becomes important at small M, and the 
cross section decays weaker than Although the parton model results are supposed to 

coincide with the predictions of the light-cone approach different approximations are used 
and the amount of disagreement between the results may serve as a measure of theoretical 
uncertainty. Nevertheless, we believe that light cone prediction is more reliable. 

The data in figs. 2 and 3 have been measured in a pseudorapidity interval of 2.1 < 77 < 
2.65, while our calculations have been performed for rapidity y = 2.4. Note that the data are 
also subject to pt cuts, which exclude very small transverse momenta, pt < 50 MeV. These 
cuts are not included in our calculation. In the region of our interest (M = 0.2-0.8 GeV) 
however both approaches underestimate the experimental data by at least a factor of 10. We 
therefore believe, that also a more careful calculation would not change our results. 
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3 Secondary interactions of radiated "prompt" gluons 



The second process, labeled "prompt gluons", considers gluons which are radiated in an 
elementary NN collision and which convert into a virtual photon via gluonic Compton scat- 
tering on another nucleon. To estimate their contribution to the dilepton pair production 
we start with the elementary Compton subprocess G q ^ ■y* q (see the upper part of Fig. 03 
where the on-mass-shell gluon is marked by a dashed line). The invariant mass distribu- 
tion for from this process can be expressed in terms of standard kinematic variables 
s = {pi + PgY, t = iPi ~ ViY ^-iid u = {pg — p-yY follows (we use expressions from 



with substitutions s =^ s — mi, i =^ i—mj,, u =^ u — m^] 



9' 



£a{qG^in-X) 
dM^di 



9M2(s 



m 



2\2 



m 



2\2 



+ {i-mlf + 2M 



m,. 



m?){t 



(17) 



If one neglects the transverse momenta of the initial particles, the cross section for the lepton 
pair production by the gluonic Compton process on a quark takes the following form when 
going from the variable t to rapidity y (for notations see Fig.[^) 



d^a{qG^in-X) 
dMdy 



2M 



di{y) 



dy 



d^ajqG^l+l-X) 
dM'^di 



where 



i{y) 



is + M- mi 



El — "Pi tanh(?/) 



Ei^Eg- {A + Pg) tanh(?/) 



(19) 



To obtain the cross section for a nucleon one has to integrate expression (|T8|) over the 
quark distributions /g and mean number of gluons n^, which are radiated by quarks from 
a projectile nucleon in an interaction with a target nucleon. To calculate the prompt gluon 
spectrum we use the perturbative evaluation for the cross section of gluon radiation [|1^ 



dV(g^gG) _ 3«,C2m2a4A;2 + [l + (l-a)2] (A;^ + a^m^) 



dadk'^ 



n 



a + 



91 -a 



a 



(20) 



where k'^ = pi^, C is the factor of the dipole approximation for the cross section of a qq 
pair with a nucleon (in this energy range C ~ 3 [^) and a is the fraction of the quark light 
cone momentum carried by the gluon 



a = 



Pg \/PG + k'^+PG 



P2 



P2 



(21) 



Here light cone variables pf = Ei+ pj are used. By integrating (pO[) over k'^ one obtains the 
gluon distribution over longitudinal momentum 



dn. 



dp 



G 



dk"^ 



da 



dp, 



G 



d'^a{q^qG) 
dadk"^ 



(22) 



where 
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max 



(23) 

e=p+/Az (24) 



and A ^ 250 MeV is the pQCD parameter. Expression ( ^4]) is obtained from the restriction 
that the gluon formation length Iq should be shorter than the mean free path of a quark 
Az ^ 0.6 fm 

a(l - a)pt 



G 



< Az. 



(25) 



For each fixed quark momentum f^l eq. (|22|) gives the gluon distribution as a function of Pq 
in the interval from zero up to 

PGmax{P2) = T • (26) 



To obtain the prompt gluon distribution radiated by a nucleon with momentum p2 as a 
function of xq = Pg/P2 one has to integrate expression (p2|) with the quark distribution 
functions 



dn, 



G 



dx 



G 



/ dX2fg{X2) / 
q Jo Jo 

P2 Y. f /, 



VGmax{x2P2) dUQ ( Pq^ 

dpG -, — 5 Ug 

dPG \ P2 , 



X2 



dna 



dp, 



G 



PG=XGP2 



where 



max 



0, 



2A 



A = + Jxl + Ay pi 



These expressions give the distributions for prompt gluons in the interval < xg < xg 



2fi 



/i = 1 + Jl+mypl 



(27) 
(28) 

(29) 

•nax 

(30) 



Finally, we combine the expression ([T8|) for the cross section for the elementary process 
with the distribution function fq{xi) for the quark and the distribution dnc/dxc of the gluon 
to the cross section for the prompt ( "P" ) process on the nucleon 



'2 „P 



+ 



where 



d^a^ 
dMdy 

d'^a'' 



dMdy dMdy dMdy 

dnc d^a{qG^in-X) 



dMdy dMdy 



> / dxifq{xi) / dxG 
„ Jo Jo 

{y -y) ■ 



g 

d^a^ 



dx 



G 



dMdy 



(31) 

(32) 
(33) 
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For nucleus-nucleus collisions the geometric factor has to take into account that points 
of gluon creation and interaction are different 



AB 



dMdy 



/r roo roo 

dh Ids / dzAPAi^.ZA) / dzB Pb{^-^,zb] 



(34) 



X a 



NN 



dzpA{s,z) ~ + / dzpB(h — s,z] 



dMdy 



2s 



dMdy 



Again, expression (0) has to be divided by the total inelastic cross section to compare 
with experiment. Our results for proton (p-Be, p-Au) and heavy ion (Pb-Au) scattering are 
shown in Figs. ^ ^ 



4 Conclusion 

We have calculated two contributions to the spectrum of invariant masses of dileptons pro- 
duced in proton-nucleus and nucleus-nucleus collisions. The direct production of dileptons 
NN — » l^l^X and the production by GN l~^l~X from gluons which have been produced 
in another NN collision. Both processes happen during the early phase of hadronization and 
may be treated by pQCD on the partonic level - in contrast to the final stage of the nuclear 
collision where hadrons and their decay into dominate. 

The results of our calculations are summarized in Figs. |^ and |^: the mechanisms consid- 
ered contribute to the observed spectrum less than 10% and therefore are unimportant on 
the present level of discussion. On this level also the theoretical uncertainties, for instance in 
the choice of the value for constituent quark mass or in the evaluation via two mechanisms 
(light cone vs. parton model) are not yet of importance. 
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